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Abstract 



A full analytic calculation of the two-loop diagrams contributing to the static 
Q^l potential in QCD is presented in detail. Using a renormalization group improve- 

Q^' ment, the "three-loop" potential in momentum space is thus derived and the 

(D . third coefficient of the /3-function for the ^-scheme is given. The Fourier trans- 

formation to position space is then performed, and the result is briefly discussed. 



X. 

c^ ■ 1 Introduction 

The famous Coulomb potential of electrodynamics is very important as an essential 
ingredient in any non-relativistic problem involving charged particles and consequently 
as being responsible for many phenomena in everyday life. It is thus no surprise 
that its analogue in chromodynamics has also been of great interest since 20 years. 
Although the QCD potential certainly is not as ubiquitous in the macroscopic world, 
it still represents a fundamental concept which, besides the fact that potential models 
have been astonishingly successful in the description of quarkonia, might give us some 
deeper understanding of non-abelian theories and especially of confinement. There is. 
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of course, no known way to analytically derive the confining part of the potential from 
first principles up to now, and consequently this paper will also be restricted to an 
analysis of the perturbative part. Nevertheless one may hope to obtain some hints on 
the non-perturbative regime in this way. In addition, a calculation of the perturbative 
potential is required for comparisons between continuum QCD and lattice results. 

The first investigations of the interaction energy between an infinitely heavy quark- 
antiquark pair date back to 1977 [jl], |^, H and were initiated by L. Susskind's Les 
Houches lectures [Q. In these works some general properties of the static potential were 
discussed and the one-loop and parts of the higher-order diagrams were computed. In 
the years that followed, much effort went into deriving (quark-mass suppressed) spin- 
and velocity-dependent corrections in various frameworks, but it is surprising that 
for a very long time there was no complete two-loop calculation of the true static 
case available. Only recently has this gap been closed [Q, and the present paper is 
essentially an extended version of [Q presenting details about the two-loop calculation, 
which might be of interest for other problems as well. 

Before turning to the real problem, however, a brief description of the structure of the 
paper should not be missing: as an introduction of the notation and for illustrative 
purposes, the QED potential is discussed in section ^, which is followed by a section 
describing the main additional problems encountered in the non-abelian theory. Section 
H is concerned with the techniques needed for the actual two-loop calculation, and 
finally sections ^ and ^ contain the result for the potential in momentum and position 
space, respectively. 



2 Abelian case: the QED potential 



A definition of the static potential in QED which is both convenient for analytical and 
lattice calculations and which in addition is manifestly gauge invariant can be given 
as follows: 

consider the vacuum expectation value of a rectangular 
Wilson loop of spatial extent r and temporal extent t, 



(0|T exp{ieidXf,A>'}\0). 




If we let t approach infinity, the lines corresponding to |t| = oo will only give a negligible 
contribution to the line integral because the length of the integration path is much 



smaller than that of the other two line^. In other words, 
(0|T expiieidx^Af^JlO) 
'^ (0|T exp{te f'^dr{AoiT,-r/2) - AoiT,r/2))}\0). 



The path integral representation of the remaining expectation value equals the partition 
function of a system described by the usual (purely photonic) QED Lagrangian with 
the addition of a source term J^{x)A'^{x) with 



Jf,{x) = ev^ 



*'»)(x+i)-.<»>(x-i) 



Vfi = S^o- 



As the partition function is dominated by the ground state energy in this limit, i.e. 
the path integral approaches exp{—itEQ{r)), and as the ground state energy is exactly 
what we would term the potential, we are led to define 

V(r) = - lim -ln(0|T expiie /rfx.A^llO). (1) 

t^oo it J 

For pure QED the vacuum expectation value is a gaussian path integral and an exact 
calculation is therefore possible, with the result 

where as the second term the expected Coulomb potential appears, but there is also 
an infinite constant representing the self-energy of the sources known from classical 
electrodynamics. 

An exact solution will no longer be possible as soon as light fermions are included or 
a non-abelian theory is considered. It is therefore useful to perform a perturbative 
analysis as well. The Feynman rules for the sources can be read off easily: each source- 
photon vertex corresponds to a factor iev'^, an anti-source obtains an additional minus 
sign. When we expand the time-ordered exponential, we may introduce 0-function 
to express the different possible time-orderings of the fields A'^, and in turn can re- 
interpret them as source propagators in position space, 

Sf{x - x') = -ie{xo - a;o)5(^)(x - x'), (3) 

the reverse time ordering must be used for the anti-source. Transforming the expression 
to momentum space one obtains 

Sf{p) = ^- (4) 

vp + le 

^It should be evident from the definition (|l|) that a small contribution from these two lines would 
indeed not affect the potential. 




a) 



Figure 1: One-loop ladder diagrams. The double-lines represent the static sources. 

with V -^ —V for the anti-source. One should note that in the path-integral approach 
we could in fact omit the time-ordering prescription, as it is already implicit in that 
formalism, but the final integrals we would encounter would not be easier to solve. 
The approach we take, however, will prove to be very useful in the non-abelian theory. 

In addition, an immediate observation is that the Feynman rules are exactly the same as 
those of Heavy Electron Effective Theory, the QED analogue of HQET, where v would 
represent the electron's velocity. This comes as no surprise, as we are investigating the 
infinite mass limit of QED, and it is well known that the potential can also be derived 
from the scattering operator. The static QED potential thus should be derivable from 
the scattering matrix of HEET, correspondingly the QCD potential from HQET, and 
this is in fact what is done in practice. This approach can even be used to determine 
the spin-dependent corrections to the potential 0. 



At tree level, one obtains the scattering amplitude 



ie dxodyo v^VyD^''{x - y) 



2 / d^q e^'i'" 



ieH 



x=0,y=r J (27r)3-q: 



i2 



where D^^ is the photon propagator. At first glance this seems to be the complete 
result already and one might wonder what happens with the higher order diagrams. 
But some care is required because the sum of all diagrams does not correspond to the 
potential, but to its exponential. That is exactly what the loop graphs are needed for. 
At one-loop order for example, the ladder diagrams shown in Fig. |I] appear. When 
working in position space it is easy to see that summing the two effectively removes 
the anti-source propagator — because it is only a G-function. Hence, adding them 
once more with x ^-^ x' the source propagator is also removed and one finds 



2! X Fig. I = ( - ieH j dxodyo D°°(a; - y) 



x=0,y=r 



(A similar analysis in momentum space is also feasible, but more difficult [Q.) It is 
clear from the way the ladder diagrams are constructed that this behaviour persists to 
all orders: the ladder diagrams are derived from the uncrossed ladder by permuting all 
vertices, i.e. generating all time-orderings, on one of the two fermion lines. Hence after 



summing all of them we end up with the uncrossed ladder again, where the G-functions 
on one of the fermion lines are removed. The source propagators on the other line are 
removed by adding the same sums, with all possible permutations of the names of the 
vertices, just as was done in the one-loop example. Since there are n\ permutations at 
n-loops, the exponential thus forms. 

The other types of graphs — source self-energy and vertex correction diagrams — 
can be shown, along the same line, to produce the products of the form E" • V^^^ (r) 
predicted by the exact formula (§). 

If light fermions are to be included, an explicit exact expression for the potential can 
no longer be derived and one has to rely on the perturbative treatment. The effect 
of the fermions is to introduce a running coupling in two different ways. The obvious 
point is that fermion loops cause vacuum polarization effects, which can be accounted 
for by defining an effective charge in momentum space via 

ttcff(q ) 



1 + n(q2) 

where n(q^) represents the vacuum polarization function and a the fine structure 
constant. This definition is gauge invariant and includes a Dyson summation. One 
might thus guess that 

q 

but the formula need only be correct at the one- and two-loop level, because starting at 
three loops light-by-light scattering diagrams enter. The correct formula should read 

^ = _±L^, (6) 

q^ 

with Q^v 7^ Ooff, which, however, merely serves as a definition of a^. Nevertheless, this 
definition is quite convenient, especially when turning to the non-abelian case where a 
gauge-invariant definition of an effective charge is nontrivial. 



3 Complications in the non-abelian case 

One novelty that arises in QCD is well known: the non-linear nature of non-abelian 
theories introduces additional diagrams which cause a running coupling even in the ab- 
sence of light fermions, and the fact that gluons carry colour also requires a redefinition 
of the Wilson loop and consequently of the potential as follows: 

Y[t) = - hm -ln(0|Tr Pexp (ig jdx^ A^T")|0). (7) 



The generators T° have to be inserted in order to absorb the gluons' colour index, 
and the fact that the generators do not commute requires the introduction of a path- 
ordering prescription. An important point is that this prescription is not automati- 
cally taken into account by using the path-integral formalism, in contrast to the time- 
ordering. Consequently, all we can do even in the path-integral approach is use the 
definition of the path-ordered exponential, i.e. expand it, and the usefulness of the way 
we analyzed the abelian case should become clearer. The net effect of the modifications 
in the non-abelian case is a complication of the way the exponentiation of the potential 
arises, and that the sources are automatically in a colour singlet state. 

In principle we are free to choose any representation for the generators, but since we 
intend to describe quarks we are going to use the fundamental one. On tree level, the 
only change with respect to QED then is the colour factor Cp = Tp^N"^ — 1)/N that 
multiplies the coupling in the potential, where A^ is the number of colours and Tp the 
normalization of the generators. It is thus convenient to define, in analogy to (|^), 

47rav(q^) .„. 



Viq') = -C, 



q2 



as this has the following advantage: Cp is the value of the Casimir operator T^-T"- = 
CrI in the fundamental representation T°- = A"/2, where the A" are the Gell-Mann 
matrices. If we chose the adjoint representation T^j = — i/"*-', which would be equally 
acceptable, the only thing that would change in (^ is the overall factor Cp which 
would be replaced by Ca, the coupling cty would remain the same. The statement is 
trivial at tree level, but it is in fact true to all orders as will be explained after the 
discussion on the exponentiation. By choosing the adjoint representation we obtain 
the static potential for gluinos in a colour singlet state. 

In higher orders the presence of the generators causes the individual diagrams to obtain 
different colour factors and consequently the exponentiation has to be more involved 
than in the abelian case. Although this point has already been analyzed in [Q , it seems 
appropriate to recall the discussion here. It should, however, suffice to consider the 
ladder diagrams to explain the basic idea. 

The QCD one-loop ladder diagrams corresponding to Fig. |l| obtain the colour factors 

Fig. la oc C| , Fig. |l|b oc C| - ^Cp (9) 

and we immediately can identify the abelian-like terms oc Cp that are needed to build 
the exponential of the Coulomb potential. But it is also obvious that the remainder of 
Fig. |I|b, together with a contribution from the vertex corrections involving the same 
colour factor, gives an additional contribution to the potential, which then has to be 
exponentiated by the higher order diagrams. 

Fig. 1^ shows the two-loop ladder diagrams, which involve the colour factors 

Fig. 2a oc Cf, (10) 





a) 



b) 



c) 






d) 



e) 



Figure 2: Two-loop ladder diagrams in QCD. 



Fig. |6,c oc C|- 


~^^ 2 


Fig. Id, e oc Cp- 


- 2C|^ + ^Cp 


Fig.0/ oc C^- 


- 3^1^ + ^Cp. 



(11) 

(12) 
(13) 



As expected, each diagram contains the "abelian" term C%, thus forming an iteration 
of the tree-level potential, but the terms oc C%Ca are also iterations, exactly those that 
arise from exponentiating the additional one-loop term. This can be seen as follows: 
when working in coordinate space and neglecting the colour factors, we can look for 
combinations of the diagrams in which one of the anti-source propagators is removed, 
with the result 



ih) 



(/) 



(c) + {d) + (/) 



{d) + 



if) 




Summing the three lines we find that the coefficients of —C'^Ca/'^ in Eqs. (|TT]-0) 
equal the corresponding diagrams' coefficients in the sum and that the B-functions 
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fixing the position of tlie marked gluon line are removed. Hence tlie sum is a product 
of one-gluon exchange and the crossed ladder diagram with exactly the right colour 
factor. A similar analysis for the remaining two-loop diagrams leads to the conclusion 
that the exponentiation works in the non-abelian case as well, and that all terms oc Cp 
or CpCA indeed constitute iterations. 

In consequence, the number of diagrams that really have to be computed to determine 
the two-loop contribution to the QCD potential is somewhat reduced. Nevertheless, 
when compared to QED the number of diagrams is still quite large, as in the abelian 
case only very few remain (essentially the vacuum polarization). 

Having discussed the exponentiation we are able to proof the statement that ay is 
independent of the representation the sources are in. It is obvious that for any repre- 
sentation R, the colour factor for a given diagram is the same as the one for R = F 
if we replace Cp by Cr, because it is completely determined by the algebra of the 
generators: only C/j and Ca can appear since the generators themselves and the struc- 
ture constants arising from commutators are the only "matrices" around, and their 
combination is fixed. We have also seen that at each loop order the only net contri- 
bution to the potential arises from the colour factors linear in Cp, all other terms are 
iterations. The definition (|^) then implies that a^ does not involve Cp at all^ and thus 
the replacement Cp -^ Cn does not affect the coupling. 



4 Details of the two-loop calculation 

It is worth presenting at least some details of the techniques used to compute the 
two-loop diagrams, first of all because they are helpful when trying to reproduce the 
results, but also because they might prove useful for other calculations as well. 

The remaining diagrams are best analyzed in momentum space, using the kinematics 
that follows from the Wilson loop definition: the "on-shell condition" for the source 
reads vp = where p denotes the four-momentum carried by the source. Thus the 
sources may have any three- moment a, the actual values of which are, however, unim- 
portant, as the only quantity that appears is the momentum transfer q^ = (0, q). 

By choosing a convenient gauge the number of graphs can be further reduced: if we 
employ Feynman gauge, all diagrams with a three- or four gluon vertex with all gluons 
directly coupled to source lines vanish. This welcome feature is caused by the replace- 
ment of the Dirac matrices 7^ in the source-gluon vertex by the simple factor v^ = 6'^^, 
which means that all vertices are interchangeable as far as their Lorentz structure and 
momentum dependence is concerned (the momentum dependence is mentioned because 

■^We only consider factors Cp coming from the sources here. Of course a^ does involve Cp, but 
these terms arise from fermion loops and are unrelated to the representation of the external sources. 
Consequently they should not be replaced. 



it is the property which is destroyed by other gauges). The symmetry properties of 
both the three- and four gluon vertices then imply that the diagrams vanish. 



rs^'^^^^l^^^^^ r^TS'iS'^^^^^?^^ rf^^^^^^^ 





a) 



b) 



c) 






d) 



e) 



Figure 3: Two- gluon exchange diagrams in QCD that remain in Feynman gauge. 

Consequently, apart from the ladders Fig. 0d-f, the two-loop gluon self-energy and 
double insertions of one-loop corrections, only a few of the two-loop vertex correction 
graphs and the six two-gluon exchange diagrams of Fig. ^have to be calculated. In fact 
the number of diagrams is even slightly lower because Fig. ^ has a vanishing colour 
factor for sources in the colour singlet state, and of the ladder diagrams only Fig. ^ is 
required, because the relation 



Fig. 2d+|e = -Fig. 



(14) 



holds if the corresponding colour factors are neglected. The equality can easily be 
proven in momentum space by applying the trivial identity 



1 



1 



1 



1 



1 



1 



Iv — kv + is kv + ie 
to one of the two source lines. 



Iv — kv + ie Iv + ie kv + ie Iv 



le 



Using dimensional regularization with D = 4 — 2e for both the infrared and ultraviolet 
divergencies, the integration by parts technique can be used to reduce most of the 
integrals that arise to products or convolutions of four types of one-loop integrals, that 
in turn can be computed by standard methods. With the notation 



~_ (27r^^ 



vn^ 



Vn{v) = 'f^v^'' ■■■T^ 



G{n,m;k,j] 



T{n + m — j — ^ 



{A7c^y ^'''^''\ / ^' B[ --n + k-j,- 
1 ml m 



D 

2 



D 
2 



m + j 



(15) 



2 r(f -n + k- j)T{2n + m - k - D) 



GH{n,m;k,j) = (4^)^^^ ^ \;, , '- (16) 

i [n)L [m) 

the integrals are, omitting the ie in the propagators for brevity: 



[fc/2] 



\—p J \pv + w J 



(17) 



dpi 



1\"/ 1 \"^/1 



'P^ J \~{p ~ Qy J \~p'v , 



iig=0 



9x£_^_„_m VVT 



rfp 



V)^— "^;^G(n, m; -a, -a/2), (19) 

-,\n/ \ a, / \ b 

I \ I W \ I W \ 



'P J \pv + w J \pv J 

-2wr-'- ^^^~\~^^ Gn{n^ a- -b, -h). (20) 

i[D — In) 



The formula for the standard massless two-point function Eq. (|T7|) is a well known result 
which can be found for example in ^ . Its HQET analogue (0) and the HQET three- 
point functions ( |19|) and (pO]) can be calculated in a similar way using the parametriza- 

tion 

1 _ T{n + m) /-^ a^-^ 

'cfHf^ ~ T{n)T{m) Jo " (a -F «&)"+'" 

to combine gluon and source propagators, instead of the usual Feynman parameters. 
Eq. (^) without the factor Vk{p) has already been given in P], and Eq. (|20|) can be 
derived from Eq. (8) of 0. 

The reduction of a given graph to the two functions G and Gh can be automated 
with the help of a computer program such as FORM [^, but an additional final 
reduction to a few basic (j(j:/)-functions seems quite difficult because of the appearance 
of e-dependent values for the last two arguments of G and Gh and additional ratios of 
F-functions. As an example consider the integral 

dpdp' 



p'^ip + qY{p'Y{p'vY{p'v + pv){pv) 
^ 2^-2F(D-4)^ .. . o ON /-X 1 



-GHil,l;-2,-2)Jdp- 
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m 



-8r(-2e)r(l + e) 



Gh{1, 1; -2, -2)^(1, 1; -2 - 2e, -1 - e){-q' 



-l-2e 



r(2-2e)r(2 + 2e) 

which corresponds to Fig. ^. It thus seems easier to directly rewrite the resulting 
expressions in terms of F-functions and immediately expand them in e. 

However, as already mentioned, not all of the graphs can be reduced in an easy way 
to the above one-loop integrals, a few diagrams remain that involve the following type 
of true two-loop integrals: 

I(a,b,c-n,m)= flplr(^] ( ^] ( j-^L_] ^^ (21) 



"5 ^^5 ^5 '"5 



p. 



r"^ ) \{p ~ r — qY ) (pv)"' {rvy 



with vq = and n,m > 0. From 



d 



dpdrVf,-—f{p,r;q) = 
dpf, 

one can derive the recursion relation 

nN+J = 2(aA+ + cC+)N-/ - 2cC+M-J. (22) 

Here the operator N"*" means that the argument n oi I should be increased by one, and 
correspondingly for the other operators. The relation can be used to reduce integrals 
with n > 1 or r/i > 1 to those with n = m = 1. For example, with ra = 1, m = 2 we 
find for the integral corresponding to Fig. ^: 

/(1, 1, 1; 2, 2) = 2/(2, 1, 1; 0, 2) + 2/(1, 1, 2; 0, 2) - 2/(1, 1, 2; 1, 1). (23) 

The first two terms can be computed with the help of (p^-(pOD, but the last one 
cannot be simplified any further by the recursion relation (due to the absence of the 
propagators l/(p — q)"^ and l/(r — g)^ there are no other simple relations). 

It turns out that the following three "irreducible" integrals, the calculation of which is 
sketched in the appendix, are needed: 

11 1 

(24) 



/(1, 1,1; 1,1) = 


= |(-g^)-"G(l 


,l;-l,-l)G{ 


1, 


1 


-'-^' 


/(1, 1,2; 1,1) = 


2 Avr/i^y^" 


_ + 4 _ ^2 
e^ e 6 






e(8 ^vr2 + ^3^C3)+0(e^)] 






/(2, 1,2; 1,1) = 


1 /4vy^ 


- + - - 4 - -vr^ 
e^ e 3 






-^(24-17vr^- 


- 64C3) + 0{e 


') 







(25) 



(26) 



With these equations all the basic formulae for the determination of the two-loop dia- 
grams are given. The expressions for the individual graphs will, however, not be listed 
here, we will directly turn to the complete result instead. 
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5 The QCD potential in momentum space 

A convenient way of writing the QCD potential in momentum space is^ 

V.(q^) = -c/-^^ (27) 

Qv(q^) = an.(f.')i;s»(AVq')(^^)" (28) 



n=Q 



CxD / 2 \ 

'Ms(q)2^an(^ j (29) 



with ao = ao = 1 and where 



— Ca - -^TpUf , ai = ai + /:/o In - 
9 9 q 



2 



oi = —Ca - —Tprif , ai=ai+ po\n— (30) 



is the long-known one-loop result and 

/4343 2 TT^ 22 N. /1798 56 n^^ 

-(f - 16C3)CFTj.n, + {^Tpnif (31) 

2 2 

52 = a2 + /52ln2^ + (/5i + 2/?oai)ln^ (32) 

q^ q^ 

is the new information added by the present analysis. Note that the last term of a2 
could have been predicted from the one-loop result, it is exactly the contribution that 
would be Dyson-summed in QED by introducing an effective coupling. The third term 
originates from the two-loop vacuum polarization and thus would also be included in 
the effective coupling in the case of QED. 

Knowledge of 02 now allows us to consistently use the three-loop expression for the 
running coupling in the MS-scheme. An equation of the form (l29|) can then be used to 
derive the /5-function of one of the two couplings from the knowledge of the /5-function 
of the second scheme. If we define (3 and the coefficients /5„ in each scheme via 






-BM^-Y.l3nPf-T\ (33) 



we immediately find 



/3^(av) = l3^^{oiMs)^^-r^^ ^i*^ "ms = aMs(av) (34) 



^ The second and third equations might be incomplete for n > 2, because a„ might depend on 
In a, see the discussion in |^. The present paper, however, is restricted to the three-loop potential, 
i.e. n <2. 
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or /3o^i = /?o f and 



/?J = /?r - «i/3r + («2 - a?)/5o^s (35) 

^618 + 242(3 ll(247r2-7r^) 



(■ 
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445 + 704C3 247r2 - n* 



C 



3 



9 



+ 7, ]CiTpnf 



, 2 + 224(3 ^ ,^ .2 686-528(3 ^ ^ ^ 

+2ClTpnf + l^i_l^Cp(T^n;)l (36) 

9 

Inserting the numbers one finds that the numerical value of /3|" is considerably larger 
than that of /9^^, which, for comparison, is given by the expression |]11| 

g, ^ '^^c\ + 2ClT,n,-fcAC,T,n,-^-^ClT,nf 

+^C^{Tpnff + ^^CA{T,nff. (37) 

Hence the coupling in the potential runs faster than the MS-coupling. 

Figure ^ compares the result for a^ at the various loop-orders in graphical form, ne- 
glecting quark thresholds for simplicity: the dotted line represents the tree-level pre- 
diction, i.e. a pure Coulomb potential without a running coupling; the dashed line 
shows the one-loop result in the sense that only the (one-loop) running of Oms is taken 
into account, but the two couplings still coincide (which means that it is still tree 
level as far as the diagrams contributing to the potential are concerned and should 
thus better be termed leading order); the dashed-dotted line shows the next-to-leading 
order and the solid line the next-to-next-to-leading order results. It is evident that 
the two-loop contribution is nearly as important as the one- loop effect: the additional 
shift in Ov caused by 02 is roughly two third the size of the shift introduced by ai, and 
both corrections increase the coupling. The impact of this shift on the would-be topo- 
nium system, as an example, should be measurable: in the interval |q| = 10 . . . 30GeV 
including the NNLO-corrections amounts to a net increase of a^ by about 5 to 7%. 
The charmonium and bottomonium systems on the other hand are mainly sensitive to 
momenta below 2GeV and thus lie predominantly outside of the perturbative domain. 

At this point we should mention that for simplicity we have chosen a fixed number of 
five light flavours to be valid over the whole range of momenta. For a more realistic 
study, the bottom and charm quark thresholds would, of course, have to be taken into 
account by matching effective theories with decreasing Uf in order to restore decoupling. 
Such a procedure would, however, leave our results qualitatively unchanged, it would 
mainly cause an even faster running of both couplings at low momenta. 
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q[GeV] 



Figure 4: Comparison of the results for av(q^) at different loop-orders, where av(q^) 
is determined from Q;Ms(q^). Input parameters are aMs(M|) = 0.118, Uf = 5. 

The plot just discussed was generated by employing ( [25| ) as it stands, i.e. by evolving 
the MS-coupling, taken to equal 0.118 at the Z-mass, to the required scale (assuming 
Uf = 5) and then calculating a^ at that scale. In principle the evolution can be 
performed in different ways: one may express the coupling in terms of the QCD scale 
parameter Aqcd, or one may employ the QED-like formula 



a(g2 






'a{iJ. 



1^' 

2 ■a(/i2) 



-(^)>^ >"!.-( 



IJ- 



An 



,„lj l^„^,W') 



^^ 



/i^ 



(3J 



Although the two approaches are formally equivalent, in practice only the second one 
guarantees that the running coupling really coincides with its input value a(/i^) for 
q^ = /i^. Therefore the second approach has been adopted throughout this paper. 

One may argue that to first evolve the MS-coupling and then determine the potential is 
not the best idea because, as is evident from Fig. H, the expansion parameter becomes 
large at low scales and thus the next terms of the perturbation series might also become 
important. An alternative would be to determine a^ from (^) at the Z-mass and then 
evolve it to the scale q^. It turns out, not surprisingly, that the result is essentially the 
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Figure 5: Comparison of the running of av(q^ 
at three loops. The dotted curve displays a^ 



(dashed line) 
with the initial value av(M|) = 



(solid line) and aMs(q 



same, as can be seen from Figure |^ where the running of the two couplings is compared: 
the solid line shows the three-loop running of a^, and the dashed one a-^, again for 
aMs(M|) = 0.118 and five flavours; the dotted line shows the corresponding result 
for Q^v for the hypothetical case av(M|) = 0.118 as well, and thus really displays the 
different /3— functions. It exemplifies that the large numerical mismatch between the 
two couplings at low scales mainly arises from an amplification of the small mismatch 
at large scales. Note that in a sense the two couplings are the same at the two-loop level 
as the /9-functions then coincide and a difference merely arises because the evolutions 
start at different values. 



6 Position space 

With the momentum-space representation of Vqcd at our disposal, we are able to 
compute the real analogue of the Coulomb potential, i.e. the QCD potential in position 
space. In order to make the expressions simpler, it is convenient to introduce the 
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notation 



^(r,/x,w) = /i2"|- 



(27r)3 (q2)i+« 



(39) 



for the Fourier transform of a general power of 1/q^. Using a Schwinger parameter, 

1 1 



and a few relations between F-functions such as 



\AitXj <Xj O 



V{l + u) = ^- 



\nT(l + u)=jEU+y^^^u" for ImI < 1, 

^2 ^ 



several equivalent representations of JF can be derived, two of which will be useful in 
the following: 



J^(r,/i,M) 



(/xr)^" F(| + m)F(| - m) 
{fire 



F(1 + 2m) 
exp E 

\n=2 



47rr 



n 



(2" 



1-f-l^ 



(40) 
(41) 



where the first formula is applicable if —1 < u < 1/2, the second if 1^1 < 1/2 and 
where je denotes the Euler-Mascheroni constant. 

Adopting a strictly perturbative approach, we can use the original result leading to 
( pTl) (or in other words re-expand aMs(q^)) to find the potential in position space. 
From 



av(q^) = aMs{f?){l + 



aMs(/^^ 
Att 



■(-^'"^ 



+ ai 



+( 



Qms(/^^ 

4:71 



q 



q^ 



(/?o In ^y - {2(3oa, + f3,)\n^ + a^) + . . . 



... . . , , (42) 

where fi is the renormalization scale, we see that we need the Fourier transform of 
ln'"(/i^/q^) which we easily obtain from JF because 

^2 



In 



m f^ 



9"^ //i^ 



and thus 



d^q 



In" 



q2 q2 



^^m V q: 

/ Qm 



u=0 



(27r)'^ q^ q^ \du 
Hence the potential in position space becomes 



T{r,fi,u) 



(43) 



u=0 



V{r) 



Cf 1 1 H (^2^0 ln(Air ) + ai 



+ ( 



^)^(/5o^(41n^(;.r') + ^; 



(44) 



+2 (A + 2/3oai) ln(/xr') + 02)+... 
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Figure 6: The potential (times —r/Cp) in position space at different orders choosing 
IX = 1/r. The range in r roughly corresponds to the range in |q| displayed in the 
previous figures: r = 10"^-^fm ^ 0.002fm ~ lOOGeV and r = 10~°-^fm ^ 0.16fm 
~ fGeV. 



with r' = rexp(7£;). One should note that going beyond the strictly perturbative 
approach is quite difficult and may lead to inconsistencies as the running coupling has 
a pole and thus the Fourier transform of ( P7D does not exist. But as the pole is an 
artifact which arises because a perturbative expression is extrapolated deep into the 
non-perturbative regime, it should not be taken too seriously. 



The result (^^ may be exploited in several ways: we certainly would not use it as it 
stands because of the possibly large logarithms, but choose a scale /i that reduces the 
higher order corrections. The first and in some sense natural choice is /ii = 1/r, leading 
to the curves displayed in Fig. ^ where the effect of the loop-corrections on rV{r) is 
shown, using the same input parameters as before. As the tree-level prediction would 
correspond to a constant in this plot, it has been omitted. The graph is essentially the 
same as Figure ^ and therefore will not be discussed further. 

A second choice motivated by noticing that due to (^T]), /i will always appear in com- 
bination with r', is /i2 = 1/r'. This way we remove all terms involving 7^; form the 
coefficients. Note that these terms are not small numerically because in nth order they 
~ /?^ for all m < n as well as similar contributions from the other 



involve {2'yE(3o 
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A third choice, which is frequently encountered in other contexts as well, would be to 
select /i in such a way as to remove the first-order coefficient completely, i.e. /i3(r) = 
exp[— 7£; — ai/(2/?o)]/''", which leads to 

V,{t) = -Cp- [I + (-) [a, -a,- -a, + ^— ) j (45) 

with ^3 = aMs(/i|(r)). This approach is in effect similar to defining an effective charge 
like in QED, but one should remember that it cannot be interpreted as a Dyson sum- 
mation of the one-loop vacuum polarization and therefore there is no reason why it 
should constitute a superior choice. A real Dyson summation would lead to problems 
with gauge invariance except for some part of the fermion loop contribution. 

Yet another possibility, which is similar in spirit, would be to choose /i in such a way as 
to remove all nj-dependence from the coefficients a„, i.e. to use the BLM scale setting 
prescription |]12|. This procedure leads to 



.;|J+«i 47r +""2 ~(4^ 



^B.M(r) = -— «Ms - + ar "^:_ '' + ar' tJ.. ' I (46) 



with 



a-M ^ __c^ (47) 



BLM 



■ 133 - 396C3 247r2 - tt\ ^ 385 - 528(3 



"2 



VQ "^4/^ 19 ^^ '^ 



r^ = r"e'^' (49) 

r| ^ r" exp {^^^^ - ^^ig^f) « --e""-. ,50) 

Finally, a fifth choice follows from proceeding in analogy to momentum space and 
defining 

Vir) = -C.^^iiM. (51) 

r 

From (0) one may calculate the /3-function of the new coupling, with the result 

P^ = (3^ + y(^1 (52) 



and, after determining the initial value for a^ at Mz form (|^) (where there are again 
different choices for /i possible, we have used n = 1/r) evolve it to the right distance. 
It is evident that the appearance of •je and the C- functions makes ay differ from a^, 
and consequently the same statements holds for the /3-functions. 




r[fm] 



Figure 7: —rV{r)/Cp in the different schemes described in the text. 
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The five approaches obviously result in different predictions for the three-loop potential 
— the difference formally being of the order a^ — as is explicitly demonstrated in 
Figures |^ and |^. One should note that there is no a priori reason to prefer any choice 
of scale. The second choice /i = 1/r' might be an exception, as it consistently absorbs 
terms which only arise from the Fourier transformation and are not directly related to 
the dynamics. It suggests that the distance complementary to |q| is not r but r'. But 
still "non-dynamical" terms (n remain even in this scheme. The BLM-prescription, 
which requires more than a single scale, is also physically motivated, and the nice 
agreement between the two curves and the curve resulting from the renormalization 
group evolution (labelled "RGE" in the figures) may serve as an argument in favour 
of these approaches. 

The coupling ay is already quite large at a distance r ~ O.SGeV^^ or O.lfm, and its 
size is even strongly increased in the second and third approach due to a reduction of 
its scale implying that the non-perturbative regime starts at smaller separations. For 
five flavours for example we have fi2 ~ 0.56/r and /is ~ 0.41/r. As a consequence 
there is a large scheme- dependence of the potential in the region r > 0.07fm. From 
Figure ^ we see that the situation is even worse as the perturbation series must break 
down near r = O.lfm because the potential starts to bend down there already. We thus 
should not trust the perturbative result for distances larger than about O.OSfm, which 
is consistent with the assumption that the perturbative potential should be reliable if 
rAqcD < 0.07 ... 0.1 is satisfiedQ |1|]. 

An interesting question is whether our result gives any indication on the onset of a linear 
rise of the potential, or in other words on the existence of confinement. Unfortunately, 
the situation seems unclear in this respect. As Figure ^ proves, each new term of the 
perturbation series that is added makes the deviation from a pure Coulomb prediction 
larger and the potential more attractive. But according to the conventional definition 
(H), a linearly rising potential requires a negative coupling, the perturbative curve 
for dy thus tends into the wrong direction, and there is no sign of a nontrivial zero 
of this function. The only indication on the existence of a zero is the breakdown 
of perturbation theory at distances of the order r ^ O.lfm (or in other words the 
presence of the Landau pole) and the fact that, with increasing number of loops, this 
point appears at smaller and smaller distances. In some sense the gap between the 
perturbative and the linear part thus becomes broader. 



7 Conclusion 

We have presented results of a full NNLO-calculation of the perturbative static po- 
tential in QCD both in momentum and position space, including a description of the 
technique employed, and we have argued that the respective couphngs ay and ay are 



^For five flavours and a^^iMz) ^ 0.118 we find Aqcd ~ 210MeV. 
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Figure 8: Three-loop potential in the different schemes described in the text. 
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universal in the sense that they describe the potential for other infinitely heavy colour 
sources such as static gluinos as well. The /3-functions corresponding to the two cou- 
plings have been derived. The main result is that the two-loop contribution is nearly 
as large as the one-loop term, both make the interaction more attractive, and that the 
perturbative potential seems to be reasonably reliable up to distances tAqcd < 0.07. 
For larger values a strong scale-dependence remains and the perturbation series even 
breaks down above about O.lfm. Although the results may be interpreted as an indi- 
cation of a general tendency of higher loop corrections to strengthen the force between 
a quark-antiquark pair, to take them as a proof of confinement would be going too far. 

Starting from the results presented in this paper, it would be important to know, first, 
the impact of the two-loop contribution to the potential as a whole and, second, of 
the ambiguity due to different choices of scale on the energy levels and decay widths 
of quarkonia. The first point could be investigated in the toponium system as this is 
only sensitive to very short distances and does not suffer from the uncertainties arising 
from intermediate and large separations. The second point, however, would require the 
specification of some potential model for the latter region and thus the analysis would 
depend on additional parameters. Nevertheless, a revised study of the type performed 
in |1^ would be very interesting. 
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Appendix: Calculation of the two-loop integrals 

When trying to compute the two-loop integrals /(a, b, c; 1, 1), it is convenient to com- 
bine gluon propagators using the standard formula for Feynman parameters and to 
combine gluon and source propagators using the modified version 

1 T{n + m) f^ , a™"^ 

da- 



^m^n T{n)T{m) Jo {aa + b)''+"'' 

With this technique the two momentum integrations can be performed one after the 
other, and after some rescalings of the a-parameters one is left with 

/(a, b, c; 1, 1) = (4V)^-(-g^)^- %^,^g~,i^\ /(a, b, c) (53) 

V r(a)r(6)r(c) 

where H = a + b + c+1 and 

Pfs + i-D) ri. . roc, ^^ ,_^_ 



/(a, 6, c) = 2 — — — — — dx dy da dp x 2 (1 — x) 2 

r(L — D) Jo Jo 
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y^yb+c---l(^l_y^a-l 



_ -|D+1-E 

{a + pf + a^ - + yil-y) 

xy 



Making the change of variables a = up and /?=(! — u)p, the two integrations corre- 
sponding to a and j3 can be done, with the result 



D-3 D^_^ ^.D-2-b-c, W-X 



dx dy x^-\l - x)^-''y^-''(l - y)^-'-''-'=arctan J . (54) 

V xy 

The presence of the arctan-function makes it impossible (at least for the author) to 
compute the remaining integral exactly for all values of a, b, c and e = {A — D)/2. Hence 
only the cases really needed and only the expansion in e to the order required were 
treated. 

For /(1, 1, 1; 1, 1) this task is quite straightforward as /(1, 1, 1) is in fact finite in the 
limit e — > and its expansion in e can be computed without encountering any problems. 
But it should be noted that due to ( |55| ) / must be multiplied by 1/e to obtain / and 
thus the latter integral is divergent. It should also be mentioned that the result quoted 
in section H has not been derived in the way just described, but by constructing an 
equation involving the integral. This method will be explained at the end of the 
appendix. 

The problem with the other two integrals / is that they do not exist if we take e = 

and that it is difficult to factor out the divergence. As a first step in this direction the 

substitution 

1 

X 



1 + t^y 
can be applied to turn (E^) into 



I{i, 1,2)= / dt t-^-^' arctant ■ K,{f) (55) 

Jo 



with 



K,{f) = j'dyy-^-'%l-y)-'-^\l + eyr 
Jo 

Jo 
r(-2e)r(-l-2e) 



r(-4e) 



F(-2e,-l-2e,-l-4e;-t^ 



and where F{a, b, c; x) denotes the hypergeometric function. Although this result may 
look quite impractical, it in fact helps: one factor 1/e has been factored out, and the 
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expansion of the two hypergeomtric functions is also possible up to terms which are 
not needed anyhow. Let us first examine Ki. 



F('-2e, -2e, -4e;x) 



CO /"-i(^-_2e)2\x" 

1 

=1 



1-^E n 



n=iV,=i j'-^e ) n\ 



oo / n— 1 

1-^E no-+0(^^ 

ra=l \ j=l 



X 



X 



l-eY.- + 0{ 



n=l 



n 



l + eln(l-x) + C(e'^ 



Thus we find 



K,{e) = 1 + eln(l + e) + 0{e^)h{e) 
where the residual function ki{t^) is well behaved for t — >■ 0. A similar trick can be 



applied to K2 if we first employ the relation |]T3 



F(a, 6, c; x) = (1 — x) ^F\b, c — a,c; 
which transforms K2 into 



X 



X — 1 



K,{e) = ^^ ^5Ll ^'V + t'Y+''F( - 1 - 26, -1 - 26, -1 - 46; -^ 

1 ~r 6 



Now 



r(-46) 



F(-l - 26, -1 - 26, -1 - 46; x) 



X 



n+l 



1 + 46 



1+46 



n=l \ j= 



^1 J -46 J{n + 1)\ 



i_i±^,+.(l±^f 



X 



n+l 



+ 0{e 



1 + 46 " ■ " 1 + 46 ^^n{n + l) 
= 1 - a; + 6(1 - 46)x + 6(1 - x) ln(l - x) + 0{e^). 

By noting that -F(O) = 1 and F{a, a, c; 1) = r(c)r(c— 2a)/r^(c— a) we can even improve 
the above expansion as F must be of the form F = 1 — x + r(— 1 — 46)/r^(— 26)a; + 0(a;) 
where dix) vanishes for both x = and x = 1. Thus 



K2{t^) = {l+t 



2\2t 



r(-26)r(-i-26i 



l-6ln(l + t2 



e 



1+26 



+ 0{€')k2{t' 



r(-i-46) 

where ^2 vanishes at the origin and at infinity. 

With these approximate formulae for the Ki we can calculate the two missing integrals 
from (1^) up to terms of the order 0{e^). The divergences still present can be extracted 
by writing 



26 Jo 



/(t)-/(o)e(i-t) 



(56) 
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where the function f{t) should be regular at the origin. The new integral is finite in 
the limit e — > 0, and consequently may be expanded. One should note that in the case 
of /(2, 1, 2) there is a term t^(l + 1^)^*^ in K2 that removes the bad behaviour of the 
integrand at the origin, but introduces the same kind of divergence for large t. The 
trivial change of variables t — *■ 1/t reduces this term to the above case again. It is 
essential that we treat the whole combination as of order 1: suppose, for example, we 
would expand t^(l + t^)^^, generating 2et^ln(l +t^). Inserted into (^) and proceeding 
as described this would produce an integral of the form 

/■°° , 1,9,, arctant 
io t 

which is not of (9(1) as we would assume, but of (9(1 /e). The term also demonstrates 
that it is crucial to know that the omitted functions fcj grow at most logarithmically 
after expansion in e and thus do not invalidate our power-counting in this parameter. 

Let us now come back to /(1, 1, 1; 1, 1). Starting from the original definition given in 
section |, shifting first p ^ p + r + q, using vq = and 

1 1 



{pv + rv + ie) {rv + ie) pv + ie 



and then replacing r —* — r in the first and shifting r 
one obtains: 



rv + ie pv + rv + ie 

r — p — q in the second term. 



/(I, 1,1; 1,1) 
dpdr- 



1 



— p2)(— r^)(— (r + p + qY) {pv + rv + ie){rv + ie) 



— / dpdr 



(— p^)(— r2)(— (r — p — qY) pv + ie 



1 



rv — ie rv + ie 



Hence 



/(I, 1,1; 1,1) 
dpdr 



1 



(— p2)(— r2)(— (r — p — qY){pv + ie) 
6{rv) 



1 



rv + ie rv — ie 



— 1 I dudr 

3 J (— p2)(— r2)(— (r — p — g)2)(pv + is) 

2vr. ^ 1 f~ 6(j-v) 

In the last expression we can replace iTT6{rv) by —l/{rv + ie) (due to vq = the 
principal value part of the resulting integral vanishes) and use (^) again, or we can 
calculate the effective D — 1-dimensional integral directly. Via both routes we find the 
result given in section H. 
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